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Abstract. In our review of the blazar phenomenon, we discuss blazar models, with a 
focus on the following issues: sub-parsec jets and their environment; energy dissipation 
and particle acceleration; and radiative processes. 

INTRODUCTION 

Blazars are believed to be a sub-class of radio-loud AGNs with the relativistic 
jets pointing close to the line of sight, such that the entire electromagnetic spec- 
tra are dominated by non-thermal radiation produced in the jets. These spectra 
usually reveal two broad components, the low energy one - peaking in the IR - to 
- X-ray range, and the high energy one - peaking in the MeV - TeV range. Both 
spectral components are variable, particularly strongly in their high energy parts 
[65,19,69,61,62]. Variability time scales range from years to a fraction of a day, 
where the rapid variability often takes the form of high amplitude flares. The most 
extreme events (with the largest amplitude and shortest time scales) have been 
recorded in 7-ray bands [25,44]. Flares recorded in the GeV and optical bands in 
high luminosity blazars [68,69] and in the TeV and X-ray bands in low luminosity 
blazars [41,61,55,50] appear to be correlated. These are most likely produced by 
the high energy ends of the distributions of radiating particles, and this behavior 
suggests co-spatial production of high energy and low energy components. The 
flares are presumably produced in shocks propagating down the jet with relativis- 
tic speeds and formed by colliding inhomogeneities [52,57,9,26]. However, present 
data do not allow us to exclude other possibilities, e.g., production of flares by 
inhomogeneities flowing through and shocked in the reconfinement regions [34] or 
via collisions of jets with clouds [16,14]. 

In any case, the rapid, high amplitude variability events indicate that a significant 
fraction of the blazar radiation is produced on sub-parsec scales and this provides 
exceptional opportunity to explore jet properties in the vicinity of the central en- 
gine. In order to advance such a study, one must first understand the dominant 



radiation processes operating in the sub-parsec jets. Until about 10 years ago, the 
most viable radiation model for relativistic jets was the synchrotron-self-Compton 
(SSC) model [35,43,28]. According to this model, the smooth, polarized and vari- 
able low energy component of blazar spectra is produced by the synchrotron mecha- 
nism, while the high energy component results from Comptonization of synchrotron 
radiation by the same population of electrons which produce the synchrotron com- 
ponent. However, this rather simple picture changed after the discovery of very high 
and rapidly variable MeV-GeV fluxes in many OVV (optically violent variable) and 
HP (highly polarized) quasars by the EGRET instrument on the Compton Gamma- 
Ray Observatory (CGRO): during states of high activity, those were measured to 
exceed the synchrotron fluxes by a factor 10 or more [67,24], and it was quickly 
realized that other processes besides SSC can be more important for gamma-ray 
production [15,42,57,8]. 

Two families of new models emerged, describing the 7-ray production: these were 
ERC (external-radiation-Compton) and hadronic models. Regarding the external 
sources of seed soft photons for the inverse- Compton process in the ERC models the 
following have been suggested: direct disc radiation [15]; BELR (broad-emission- 
line-region) [57]; disc radiation scattered by the gas surrounding the jet [8]; and 
jet synchrotron radiation scattered/reprocessed back to the jet by the external gas 
[27]. In hadronic models, the 7-rays are produced via synchrotron mechanism by 
secondary electrons which are much more energetic than the electrons accelerated 
directly. They are injected following hadronic processes and accompanying pair 
cascades. The primary reaction in hadronic models involves collision of ultra- 
relativistic protons with soft photons [42] or with the ambient cold protons/ions 
[2]. The hadronic scenarios, as alternative for the SSC model, have been proposed 
also to explain TeV radiation, detected by atmospheric Cherenkov detectors in a 
number of low luminosity BL Lac objects [14,51,1,47]. 

In our review we discuss: environment of sub-parsec jets in AGNs; dissipation of 
jet energy and particle acceleration mechanisms; production of non-thermal radia- 
tion and blazar models. 

SUB-PARSEC JETS AND THEIR ENVIRONMENT 

Typical electromagnetic spectrum of luminous "non-blazar" AGNs/quasars con- 
sists of three distinctive spectral components: the UV bump and two broad (but 
less luminous) components, the IR and the X-ray one. In addition, in the optical- 
UV range, strong, broad emission lines (BEL) are present. UV bumps are pro- 
duced by central parts of the optically thick accretion disc; X-rays come from hot 
disc coronae, with possible additions from the base of a jet in radio-loud AGNs 
[70,31]; near/mid IR radiation presumably comes from dust, located in geomet- 
rically thick molecular torus and heated by radiation from the center; and BELs 
very likely originate in the disc winds, photoionized by UV radiation of an accretion 
disc [21,46,10,48]. These radiation components form dense radiative environment, 



through which jets - presumably formed in the vicinity of the central black hole 
- must pave the way out. Because of Doppler enhancement, the relativistic jets 
interact most effectively with radiation incident from the front and from side of 
the jet. At distances where most of blazar non-thermal radiation is produced 
(r ~ 10 17 — 10 18 cm), such radiation is provided by the BEL region and by the hot 
dust. Energy density of BELs on the jet axis is 
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where tbel ~ 0.31/^^,46 P c [49] and f\ is the correction function which depends on 
geometry of the BEL region. If the BEL region forms a narrow torus lying in the 
disc plane, then 

K ~ 1 + (r/r BEL y (2) 



Energy density of dust radiation is 
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where 
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and £m is the fraction of the central engine luminosity reprocessed by hot dust. 

If a luminous AGN possess a powerful jet, it is a radio-loud quasar. Observa- 
tionally, it is classified as OVV/HP quasar if viewed along the jet; as a radio-lobe 
dominated quasar if viewed at intermediate angles to the jet axis; and as a radio- 
galaxy with FR II type of large scale radio-morphology (edge brightened lobes 
with hotspots) if viewed at large angles to the jet axis [29,66]. Much less is known 
about the sub-parsec environment in very low luminosity AGNs. In particular, the 
accretion process is less-well understood; we do not know about the radiation effi- 
ciency of such objects, and to what extent optically thick plasma participates in the 
radiation processes. The low luminosity radio-loud AGNs are often accompanied 
by strong jet activity. The objects with jets pointing away from the line of sight 
are most likely FR I type radio galaxies with radio morphology characterized by 
edge darkened extended radio structures. FR I radio galaxies probably form the 
parent population for BL Lac objects [29,66], which, in similarity to the OVV/HP 
quasars, are viewed at small angles to the jet axis. BL Lac objects, together with 
the OVV/HP quasars, form the blazar category of AGNs. It should be noted, how- 
ever, that the existence of two distinct radio morphologies, FR I and FR II (and 
related bimodal classification of all radio loud objects) doesn't necessary imply bi- 
modal distribution of AGN properties. Specifically, that there is a large luminosity 



overlap between FR II and FR I objects as well as between BL Lac objects and 
OVV/HP; that there is no evidence for bimodal distribution of equivalent widths of 
emission lines in blazars; that there are some radio galaxies which on one side have 
FR I type radio structure and on another side FR II type radio structure [30] ; and, 
finally, the global spectral properties of blazars seem to form a continuous pattern 
parameterized by the total blazar luminosity [24]. 



JET ENERGY DISSIPATION AND PARTICLE 
ACCELERATION 

Rapid variability and the extent of blazar spectra to GeV and TeV energies imply 
particle acceleration in situ. The ultimate source of particle energy is the kinetic 
and intrinsic energy of a jet. The simplest and most popular model of production 
of high amplitude short term flares involves collisions of inhomogeneities moving 
with different velocities down the jet [57,9,26]. If inhomogeneities are ejected by 
the central engine with bulk Lorentz factors Ff > F s >> 1, separated by a distance 
Ar e j, they start to collide at a distance 

cAr ej 2F}F 2 S 

The collisions are followed by a formation of forward-reverse shock structures, with 
the shocked plasma enclosed between the shock fronts and moving with the bulk 
Lorentz factor of the contact discontinuity surface, which for equal rest densities of 



inhomogeneities is F ~ yFfF s (for more general cases see [13,39]). As the collision 
proceeds, relativistic particles are accelerated and produce non-thermal radiation. 
For symmetrical inhomogeneities, collisions last 

FX' FX' F 2 f + r 2 
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where A' is the radial width of the colliding inhomogeneities and v'j is the velocity 
of the faster moving inhomogeneity, both as measured in the contact surface frame. 
After that time the shocks and the injection of relativistic particles terminate. 
Thus, in case of a point source, the observer located at o b s ~ 1/T will see the flare 
which lasts 

tfi = t coU (l - (3cos9 obs ) ~ (7) 

provided electrons cool faster than the collision lasts. For the flow which is modu- 
lated such that Ar e j ~ X'/F, and Ff/F s ~few, Eqs.(5-6) give 
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i.e. during the collision the source of radiation doubles its distance from the central 
object. For sources with finite size, light-travel effects must be included, and then 
the flare time is tfi ~ (t' co ii + ^'sh/ c + a/cj/F where a is the cross-sectional radius 
of the source and X' sh is its radial width (the width of the layer of the shocked 
gas enclosed between the forward and reverse shock fronts). However, because 
\' sh ~ (2v' d / c)t' coll) where v' d is the shock down-stream velocity, and because effec- 
tively in the conical jets a < r/T (no larger area than this can contribute to the 
observed radiation due to Doppler beaming), the time scale for a source with a 
finite dimension is of the same order as for the point source. 

Efficiency of energy dissipation via collisions of symmetric inhomogeneities is 
[13,39] 
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Some portion of dissipated energy is converted to the thermal energy of the shocked 
plasma, and the remaining part is used to accelerate protons and electrons up to 
relativistic energies. Protons are quite efficiently accelerated by Fermi mechanism 
in shocks (see reviews by [18,7,38]). Time scale of the acceleration is 

(10) 

where ( is the ratio of the mean free path for Fermi scatterings to the gyro-radius, 
v' sh is the velocity of the upstream flow, and 
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is the gyro-time. The value of ( is equal to u' B /kw(k), where w{k) is the magnetic 
energy density per unit wavenumber k in the turbulent magnetic field, and u' B is the 
total magnetic energy density. For Kolmogorov turbulence spectrum, w(k) oc k~ 5 ^ 3 , 
i.e. the smallest values of ( are predicted for most energetic protons [5]. Protons 
accelerated by the Fermi mechanism in strong but non-relativistic shocks have a 
power-law energy distribution with a slope corresponding to the equal energy per 
decade of energy. 

Electrons can be accelerated by shocks as well, but first they must be preheated 
to energies j t h ~ (v' u /c)m p /m e , in order to be scattered freely over the shock front 
and energized between diverging flows, where v' u is the upstream flow velocity in 
the shock front frame. There are several processes which have been suggested 
to preheat electrons [32,40,54,6,45,17,56]. The most promising scenario is, where 
electrons are heated by waves induced in the upstream flow by protons reflected 
from the shock front. 

In general, the slope of energy distribution of preheated electrons can be different 
than the slope of accelerated protons, but those accelerated further up to 7 ^> j th 



via Fermi scatterings will have similar slope to that of the protons. The large 
difference is predicted, however, for maximum energies, which for electrons is much 
stronger limited by radiative energy losses than for protons. 



RADIATION PROCESSES 

Electron energy losses 

Radiative cooling of relativistic electrons is dominated by synchrotron mechanism 
and inverse- Compton process. The rate of radiative cooling in the source comoving 
frame, in which electrons are assumed to have isotropic momentum distribution, is 
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where u' ~ u' B + u' rad< , u' B = B' 2 /8tt is the magnetic energy density, u' rad< = 
u 'rad[ u < m e c 2 /h^} is the radiation energy density within the frequency range 
for which scatterings off electrons with energy 7 are within the Thomson limit, 
and C\ = 4o"j'/3m e c. Comparing time scale of electron radiative cooling, t' ecool ~ 
j/ldj/dt'lrad, with the collision time scale, t' coll ~ Ftfi, one can find that energy, 
above which energy distribution of electrons steepens due to radiative losses (by 
unity with respect to the injection slope) is 

7c ~ Cl Tu't fl ~ ^7(r/15)(t /i /lday) ( } 

Radiation energy density is contributed by both local emission and by external 
radiation fields. The dominant local contribution to u' rad< for 7 ^> 1 is provided by 
synchrotron radiation: 

/ L syn< (u) 

where L syn< is the apparent luminosity of synchrotron radiation within the Thom- 
son limit {y < Tm e c 2 /h'y). Contribution from the external radiation fields to u' rad 
is 



u' ext = -Ji'dn' = -J iv^dn (15) 



where / is the intensity of the incoming external radiation, T> in = l/r(l — f3 cos 9i n ), 
and 8 i n is the angle between the jet axis and an incoming ray. Provided that the 
dependence of / on 6 in is much weaker than the dependence of V in on 9 in , the 
u' ext is dominated by those external radiation fields which at a distance of the flare 
production contribute significantly from the side and the front of the jet [60]. Such 
radiation fields are provided by BEL region and by infrared radiation by hot dust, 
and from Eq. (15) 



f 2 T 2 u ext (r) (16) 
where for the external sources located in the narrow torus at r ext 

h * ( 1 - -T===) ( 17 ) 

Noting that distance of the flare production, 

r n ~ ctfiY 2 ~ 6 x 10 17 (t /i /lday)(r/15) 2 cm (18) 

is on the order of r BEL , the value of f 2 can be sometimes significantly lower than 
unity. For scatterings of IR radiation, because r dust rfi, we have f 2 — 1. 
Noting that v BEL ~ 10 15 Hz, v m ~ 3kT dust /h ~ 6 x 10 13 T 3 Hz, and that 
v' ext ~ V~^v ext = \^f2Tu ext , one can find that BELs are scattered within the Thom- 
son limit by electrons with energies 

m e c 2 7.9 x 10 3 

7 < yfJ£hv BEL ~ v^(r/15) ( } 

while the dust radiation is scattered within the Thomson limit by electrons with 
energies 

m e c 2 1.3 x 10 5 

7 < TrW ~ T 3 (r/i5) • ( } 



Proton energy losses 

r via synchrotron and 

c 27 V (21) 



Protons lose energy via synchrotron and inverse Compton radiation at a rate 

dip 



, eft' 

where c 2 = (m e /m p ) 3 c 1 . 

Comparing the time scale of proton radiative cooling, £p COo; = r y p /\d'y p /dt'\, with 
the flare time scale, one can find that only protons with energies larger than 

(rnrt 1 15 X 1011 1 / x 

ry( rad ) ~ ~ (02) 
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can radiate efficiently. 

Protons can also lose energy via inelastic collisions with photons and cold ions. 
Collisions with photons can lead to meson production (photo-meson process) and to 
direct e + e~ pair production. Time scales of proton energy losses in these processes 
are [4]: 



f(e)> 



6.7 x 10 

W 



19 



(24) 



2.2 x 10 15 , . 

t' PP * —, (25) 

for photo-meson production, pair production, and pp collisions, respectively, where 
n^' = J (tt)i n' u du f , n>' = J ( e y n' u du', and the threshold photon energies for pj 

V th V th 

collisions are v$' ~ m n c 2 /h^ p ~ 3 x 10 22 /7 P , and v^' ~ 2m e c 2 /h^ p ~ 2x 10 20 /7 P . 



BLAZAR MODELS 



Depending on the dominant mechanism considered for the 7-ray production, 
blazar models can be divided into three groups: 

- SSC models, where production of 7-rays is dominated by Comptonization of 
locally emitted synchrotron radiation; 

- ERC models, where production of 7-rays is dominated by Comptonization of 
external radiation fields; 

- hadronic models, where production of 7-rays is dominated by synchrotron radia- 
tion of proton-initiated-cascades (PIC), and/or by synchrotron radiation of protons 
and muons. 

These models are discussed below as applied to OVV/HP quasars and TeV-BL Lac 
objects. 



SSC models 

The simplest SSC model is the one-component version, where the radiation ob- 
served at any given moment is produced by a source/shock moving along a jet. In 
the SSC model, the frequency of the 7-ray luminosity peak is located at 

(7>l if hu L > lm m e c 2 /T 

V~f m m e c r/h otherwise 



where 



vl ~ c Bl 2 m TB' (27) 



is the location of the synchrotron peak, 7 m is the location of the break in the 
electron energy distribution, and cb = (2/3ir)(e/m e c). These equations together 
with the ratio of the luminosity peaks, 

v hL Vh ^ u' syn < 
v L L VL ~ u' B 



allow to derive B', T, and j m . 
SSC in OVV/HP quasars 

The SSC process can dominate 7-ray production in OVV/HP quasars if 
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This can be satisfied in OVV/HP quasars only if geometry of the BEL region is such 
that j'2 < 1CT 2 and if temperature of dust is lower than 300 K, or that T < 5. First 
condition is difficult to satisfy, because from reverberation mapping, tbel ~ f'fi\ 
second condition is in contradiction with observations of near/mid IR radiation in 
AGNs of radio-lobe dominated quasars. Finally, T < 5 contradicts with VLBI data. 
Furthermore, one-component SSC model predicts X-ray spectra which are too soft. 

SSC in TeV-BL Lac objects 

Much more promising are SSC models as applied to BL Lac objects, in particular 
to the BL Lac objects which are TeV sources [61,63,36,3]. In these objects, Klein- 
Nishina effects are important. Using Eqs. (26-28), one can find that typically B' ~ 
0.3 Gauss and 7™ ~ 3 x 10 4 [63]. Using these numbers, Eq. (13) gives 7 C ~ j m . This 
suggests that spectral breaks at peaks are due to effects of cooling; furthermore, 
the fact that 7 m is only by a factor few lower than 7 max ~ /iz/^/rm e c 2 explains 
why spectral components in TeV blazars are hard almost up to highest frequencies. 
Finally, noting that maximum electron energies presumably are determined by 
the balance of the acceleration time scale and radiative cooling time scale, having 
t'coii ~ t'cooi provides natural explanation for existence of "histeretic" variability 
patterns observed with both signs in the X-ray band in the spectral index vs. flux 
diagram [61,62,36,23]. 



ERC models 

In the ERC one-component model the high energy peak is located at 

VH ~ £ir 2 7^ e * t (30) 

while ratio of the peak luminosities is 



(31) 
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ERC in OVV/HP blazars 

When applying ERC model to Q-blazars, it is better not to rely on relation (27). 
This is not only because v L is located in these objects in the spectral range which 
observationally is poorly covered (far-IR band), but also because there are some 
observational indications that far-infrared radiation can be strongly contaminated 



by radiation produced at larger distance in a jet than the high energy flares. These 
indications are: the low energy break of the synchrotron component is at lower 
frequencies than the break due to synchrotron-self-absorption on sub-parsec scales 
[57]; and there is clear trend of decreasing amplitude of the flare with decreasing 
frequency in the synchrotron component [19,20,69]. However, in order to close 
a system of equations describing the model, one has to replace one equation by 
another. Such relationship can be derived by assuming that 7 m = 7 C [9]. Another 
useful approximation is to treat as a known observable T, instead of u ext . This 
is because due to unknown geometry of the BEL region and unknown maximum 
temperature of dust, model errors which can result from using uncertain values of 
f'2 and Tdust can be larger, than resulting from uncertain value of Y. Alternatively, 
instead of treating V as known, one can make an assumption about equipartition 
of magnetic fields with relativistic particles. Our preliminary results show that for 
T ~ 15 both approaches give similar results. Using observables typical of OVV/HP 
blazars, namely: u H = 10 21 Hz; u H L VH / u L L UL = 10; tfi = 1 day; and T = 15, and 
assuming that the high energy luminosity is dominated by ERC(BEL), we obtain 
B' ~ 1 Gauss; 7 C ~ 70; ~ 2 x 10 11 Hz (which actually is below the synchrotron- 
self-absorption frequency); "i ma x ~4 x 10 3 ^ z/ S2/njmaa; /10 15 Hz; and fi ~ 0.2. 

Model parameters derived in this manner can be used as "first approximation" 
parameters for dynamical models intended for a study of time evolution of electron 
energy distribution. Such models have been recently developed using two kinds of 
approximations: one is a spherical source approximation by [12,26] and another is 
a thin shell approximation by [9]. Results of both demonstrate the ability of ERC 
models to explain the observed spectra and flares of OVV/HP quasars. In addition, 
these models support the earlier predictions that the SSC radiation component, 
despite its lower luminosity than the ERC component, can dominate the soft/mid 
X-ray bands [33,37]. 

ERC in TeV-blazars 

Very little is known about radiative environment in these objects and, so, energy 
density of external radiation fields must be treated as a free parameter. But there 
are interesting constraints regarding the level of the IR radiation from dust. This 
radiation must be very weak, otherwise TeV 7-rays would be absorbed in the pair 
production process [11]. However, production of TeV component by Comptoniza- 
tion of radiation of dust is still feasible, provided T > 15. 

Hadronic models 

Hadronic models, which are motivated by theories of particle acceleration, suffer 
difficulties in explaining the observed electromagnetic spectra and/or short term 
variability. Proposed for the OVV/HP quasars, models involving proton induced 
cascades (PIC) [42] predict X-ray spectra which are softer than observed. They 
also require fine-tuning of model parameters in order to have the location of the 



radiation deficiency (the "dip" between the low energy and high energy peaks in 
the broad-band spectra) located at the right frequency in all OVV/HP blazars. 
Furthermore, there are already several examples of OVV/HP blazars, where X-ray 
spectra have slopes ax < 0.5 [22,53,64], which cannot be explained in terms of 
the synchrotron radiation mechanism. This is because synchrotron radiation in the 
X-ray band requires such energetic electrons, that they must cool efficiently and, 
therefore, produce spectra with ax > 0.5, even if injected mono-energetically. 

From unknown us reasons, no contribution of pp— > pe + e~ was taken into account 
in PIC models. The process is about 100 times less efficient per photon than photo- 
meson one, but this is compensated by a much larger number of photons above the 
threshold energy, which for pair production is ~ 100 lower than for photo-meson 
process. 

Another hadronic model proposed for OVV/HP blazars is based on the process 
of inelastic collisions of ultrarelativistic protons with the cold ions. Such process 
in order to be efficient requires very dense targets and huge bulk Lorentz factor 
(T ^> 10) of streaming protons. In the paper [2] it is proposed that the dense 
target is provided by the walls of the funnel formed by geometrically thick disc. 
The main problem with this model we can envisage is overproduction of soft X-rays, 
which would be produced due to Comptonization of accretion disk UV photons by 
cold electrons in a jet This effect can be avoided if in the vicinity of the black hole, 
up to at least 100 gravitational radii, jet is not yet collimated [59]. But then beamed 
blazar radiation must be produced at larger distances, probably at r ~ 10 17 ~ 18 cm 
as the flare time scales suggest. 

Hadronic models have been proposed also for TeV-blazars. In these objects the 
soft spectra of PIC models contrast with the very hard observed spectra even more. 
This problem is avoided in proposed recently the proton-synchrotron model [1,47]. 
In such a model high energy spectra are produced directly by protons, just via 
synchrotron radiation. The model predicts hard spectra and is also attractive, 
because can explain constant slope of high energy tail of TeV component during 
high amplitude flux changes [1]. In the model, maximum proton energies 7p, ma x > 

10 10 / \J fi'/lOOG are required in order to satisfy condition t' p syn {^ V}max ) < tfiF. 

Another version of the hadronic model has been suggested by Rachen & Meszaros 
[51]. In their model the bulk Lorentz factor is assumed to be rather low (r ~ 3) 
and protons with 7 Pirrtax > 10 10 lose comparable energies via photo-meson process 
as via synchrotron mechanism. Following the former, the muons are produced with 
such energies, that before decay they undergo significant synchrotron energy losses, 
dominating production of TeV peak. The weak point of both above hadronic models 
is that in order to accelerate protons up to 7 Pima a; > 10 10 energies, the quantities ( 
and v' sh /c (see Eq. [10]) must be pushed to extremes, i.e., should have values 1. 



SUMMARY 



As analysis of jet energetics and X-ray observations of blazars indicate, plasma 
in sub-parsec jets can contain from few up to tens of e + e~ pairs per proton [58,9]. 
With such a number of pairs energy flux of sub-parsec jets is still dominated by 
protons, and, therefore, the structure of shocks is determined by proton plasma. 
Then, one should expect very efficient shock acceleration of protons, with maximum 
energies as high as 7 P ~ 10 9 — 10 10 , as limited by size of the source or balanced 
by energy losses. Furthermore, very large values of ( (ratio of the mean free path 
for Fermi scatterings to gyro-radius) for electrons, 10 3 " 6 , as derived for maximum 
electron energies assuming that for them t' cool = t' acc , seem to be consistent with 
the presence of ultrarelativistic protons, provided the turbulent magnetic field has 
a Kolmogorov spectrum [5]. On the other hand, very hard X-ray spectra and short 
variability time scales put severe constraints on radiative role of protons. Does this 
contradict with theoretical predictions about very efficient acceleration of protons? 
Not necessary. 

The average proton cannot reach more random energy in the shock than E p ~ 
(r„ — l)m p c 2 , where T' u is the bulk Lorentz factor of the upstream flow as measured 
in the rest frame of the shock front. Intrinsic shocks are at most mildly relativistic, 
and therefore E p < m p c 2 is expected. Comparing this with the average energy of 
electrons which, as deduced from hard X-ray spectra of some OVV/HP quasars, 
is in the range 10 — 100 MeV, and noting that limited energetics of jets combined 
with electron emissivity implies n e /n p to be at least of the order of 10, one can 
find that fraction of energy dissipated in the shock and used to accelerate electrons 
does not have to be much lower than energy channeled to protons. This, combined 
with the fact that radiative efficiency of electrons is much larger than of protons, 
can explain negligible contribution of the proton related processes to the observed 
spectra. 
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